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ABSTRACT :In this study, vibration suppression of vehicle-bridge interaction system has been studied in 
terms of bridge dynamics. In vibration suppression, some configurations of tuned passive mass dampers have 
been considered. Motion equations of vehicle and absorbers are coupled with the motion equation of the bridge 
beam using first four mode of beam and Lagrange equations. Then the coupled equation of vehicle-bridge- 
absorbers (VBA) has been solved using the fourth order Raunge-Kutta algorithm. It is proved that the 
configuration of the absorbers that are placed on the anti-nodes of first and second mode shapes of the bridge 
beam is the most effective one of all others. 
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I. INTRODUCTION 

Vibrations of bridges under moving loads 
are vital in engineering and have been studied over a 
century. Early studies of the moving load problem 
have been carried out by structural engineers who 
accepted the moving load as a moving point force. 
Later, with the implementation of the increased 
transport speed, the subject has also been examined 
by the automotive and railway engineers. Now the 
subject, as vehicle-bridge-interaction (VBI), has 
been studied in terms of both vehicle and bridge 
dynamics with the new direction such as vibration 
suppression of the both systems. In this field, some 
of analytical studies [1,2] are valuable for the 
dynamics of the bridge beams under different 
loading cases and end support conditions. When 
higher transportation speeds are implemented, the 
dynamic amplification factor (DAF) is considered as 
important due to the inertial effects of the moving 
loads in rail and other transportation systems. 
Vibration suppression and control of the bridge 
beams have been studied by [3-6]. It has been 
reported that using Passive Tuned Mass Dampers 
(PTMD), the vibration of a high-speed train bridge 
at 300 km/h had been reduced by 21 percent. 
Another application of PTMD has been reported by 
[4] for a thin bridge under successive moving loads. 
Under a periodic excitation [10] has been studied 
vibration reduction of a thin beam by PTMD. It has 


been reported in [6] that it is impossible to suppress 
all the vibration due to the VBI, since the excitation 
frequency of the moving load is not constant. For a 
constant velocity of the load [7], for vibration 
absorbers [8] , for Timoshenko beams [9], for 

optimized damped-absorbers [10] are the other 
valuable studies in this field. Plate structures under 
moving loads are also research interests and some 
FEM of such system can be found in [11]. Resonant 
response of a beam plate of a high-speed a railway 
bridge has been optimized by [12] using passive 
viscous dampers. Vibration suppression of a plate 
structure under random excitation using PTMD and 
its Hco and H2 optimizations have been given in [13] 
and it is the early study of PTMD. It has reported 
that the most effective use of PTMD is to place it at 
a point where the mode shape is maximum at an 
anti-node. One of the other implementation field of 
the PTDM is defence science applications such as 
projectile barrel interaction in various type of heavy 
and light gun systems. For example, studies [14,15] 
have reported that the muzzle vibrations of barrels 
can be decreased at 50 percent using proper 
optimized PTMD. 

One of the biggest challenges in PTDM 
use, the influence frequency of the moving vehicle is 
not constant and when the vehicle velocity is altered, 
it is changed. Secondly, there can be more than one 
vehicle on the bridge system with different 
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characteristics. Thirdly, the moving mass case, 
taking into account the effects of mass inertia, the 
natural frequencies of the bridge beam are changed 
depending on the position of the moving load. 
Further, a VBI system is coupled and its 
linearization is difficult. 

In order to suppress vibrations of the 
vehicle and bridge beam, in this study, the idea of 
the usage of the anti-nodes for the location of the 
PTDM in [13] has been extended, and developed 
such that at least the anti-nodes of the first and 
second mode shapes may help more effectively. 
Since a general vibration, response of a vibrating 
system can be approximated by superimposing of 
the effect of each vibration mode. Using this idea, in 
this study the vibration response of the beam has 
been approximated including the effect of the first 
four modes. Then tuned mass dampers are 
considered for the anti-node of the first mode at the 
mid-point of the bridge beam, and the anti-nodes of 
the second mode that are the 0.25 and 0.75 of the 
bridge length L. The tuning of the absorber at the 
mid-point is carried out using the fundamental 
frequency while the tunings of the absorbers at the 
0.25 and 0.75 locations, the second mode frequency 
is used. Except this, the effects of one undamped 
absorber and damped absorber are also studied as a 
comparison. Moreover, the effect of the other 
position of the PTDM on beam has been analysed 
widely. From this perspective, this study can help 
both structural and automotive engineers in order to 
extend the service life of the bridges and predict the 
dynamic forces applied to the vehicle from VBI, and 
to achieve desired ride-safety and passenger 
comfort. 


II. Mathematical modelling 

In order to reduce vibrations of the vehicle 
and bridge due to vehicle-bridge interaction, a 
special system shown in Fig. 1 has been studied. 
The system consists of a single degree of freedom 
vehicle, a simply supported Euler- Bernoulli beam 
on which the vehicle moves with constant velocity v, 
and pasive vibration absorbers that are placed at 
0.25, 0.5 and 0.75 of the length of the beam. Where 
m is the mass of the vehicle, k and c spring constant 
and the damping coefficient of the suspension 
system, while y is the vertical displacement of the 
vehicle body; and parameters k, and m,(\ =1,2, 3), 
respectively, are the stiffness and mass of the 
absorbers at the given locations i= 1, 2, and 3, as 
shown in Fig. 1. Symbol y^is the displacement of 
the wheel at the contact point of the beam and wfx, 
t) is the displacement of the beam at point x and time 
t. The displacements of the absorbers are 
represented by y /, y 2 , yjfrom left to rigth. 

In the formulation for the vehicle-bridge- 
absorbers (VBA) analysis following assumptions 
will be adopted: 

• The bridge is modelled as a simple 
supported beam based on Euler-Bernoulli 
theory. 

• The vehicle is modelled as a single DOF 
lumped parameter system. 

• Only one vehicle is accepted moving on the 
bridge with constant velocity v. 

• The wheel is always in contact with the 
bridge. 

• The effect of road roughness upon vehicle 
and bridge dynamic is not considered 
during analysis in this study. 
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Figure 1. Model of Euler-Bernoulli Bridge beam with attached PTMDs and subjected to moving vehicle. 


With all these assumptions, for the system 
of vehicle-bridge- absorbers (VBA) shown in Fig. 
1 the kinetic and potential energy are expressed as 
follows, respectively: 


1 J//[vv ( , 2 (.r,0]'A + »b’ 2 (0 + 
2 ° 

mjltf) + m 2 yl(t) + m 3 y\(f) 


k ( 1 a ) 


' L 

} El [<(x,f)]& + k [y(() - w„ (f (f),r)f + 
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E.= 


K hi (0 - w»(L / 4,r)f + i, [y 2 (I) -w b (LI 2,of 
+fc 3 [3> 3 (t)-H i (3L/4,0f 
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(lb) 


W b {x,t) = Y J (pjx)T] bi {t\ 
1 = 1 
n 

w b {x,t) = Y J <Pjx)t] b jt\ 

i =1 


wf(x,f) = '£ t <p’Xx)rj bi (f), 


<Pi(x) 



i = 1,2,...,«. 


( 2 ) 

(la) 


where;//,, is the i-th time dependent 
generalized nodal coordinate,^, is i-th mode shape 
function. The applied static axle load of the 
vehicle at the contact point of the wheel can be 
determined as: 


where /2 and El, refer to the mass of the 
unit length and flexural rigitidy of the bridge 
girder, respectively. 

For the system shown in Fig.l, in order to 
obtain the equations of motion one can use the 
virtual work principle, Hamilton's principle and 
D'Alembert's principle. In this study, motion 
equation of the VBA integrated system is optained 
using Langrange’s equations and mode expansion 
method. For any point x on the beam at time /the 
deflection function w b (x, t) can be aproximated 
using the Galerkin’s method: 


f c ix,t) = -(mgH(x - ^(t))), (3) 

wherc/f(v-c(?j)is the Heaviside function. 
For the vehicle- bridge system the Rayleigh’s 
dissipation function can be expressed as below: 

D =\{ c «f>? (*’o+4 j(f) - %mMmx- m }, (4) 

In Eq. (4) c eq is the equivalent damping 
coefficient of bridge girder. In addition, for the 
given system, the Fagrangian (L = E k - Ep) is equal 
to the difference of the kinetic and potential 
energies. In such a case, for the five independent 
coordinates, the Fagrange equations can be as 
follows: 


www.ijmret.org 


ISSN: 2456-5628 


Page 3 































International Journal of Modern Research in Engineering and Technology (IJMRET) 
www.ijmret.org Volume 2 Issue 7 II November 2017 


d f dL 

dtydpjt) 


dL 


dR 


dp k (0 op k (t) 


= 0 ,k = \ 


(5a) 
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df l drjjt) 


dL dD 


drift) drift) 


= Qi , i = 1,2,3,4, (5b) 


With the state variables that are: 


p(t)={y(t)} T , ?7(t) = {/7i(t) ^(t) nff) 


( 6 ) 


And the corresponding generalized forces 


Q : = J (p t (x)f.(x,t)dx, ( = 1,2.3.4, (7) 

0 

Using the orthogonality conditions and 
Galerkin approach for displacements of the beam, 
(Eq. (2)) the equation of motion is obtained for 
VBA interaction model given in Fig. 1. In this 
case, the equation for the acceleration of the 
vehicle body and absorbers are written as follows, 
respectively: 


y = {-c[y(0 - w b (x,t )] - k\y(t) - w b (x,t)])—, (8a) 

ra 

yi={-Kiy(f)-w b tL/A,ty \}—, (8b) 

raj 

y 2 ={-^[>’(0-Wj,(L/2,t)]}—, (8c) 
ra 2 

y 3 = {-k 3 [y(t) - wJ3L/4,t)]}—, (8d) 

IUj 

Dynamic equations of the bridge as the n 
second order ordinary differential equations can be 
expressed as follows: 


N i rj i (f) + c eq ipf (x)f] i (t) + n .77, (f) 

[ f + - y(t)]\ 


+A <Pi(Zi(t))\ 


+k l [w b (L/4,t) - yjt)\ + k 2 [w b (L/2,t) - y 2 (t)] 
+k 3 [w„ (3L/4, t) - y 3 (/)] = 0 i = 1,2,3,4. 


wherezl is: 


(9) 


{1. for 0 £ t < t. 

L = I 1 (10) 

to elsewhere. 

Where fyis the time when the vehicle 
leaves the bridge. The coefficient Ais time 
dependent and is used with Eigen function 
<Pi(C(t)) in order to determine motion equation. 

In Eqs. (8a), (8b-c) and (9), shows the 


total eight second order differential equations for 
the vehicle, bridge and absorbers. The first four 
equations are for vehicle and the absorbers and 
they are transferred into first order equations by 
[16]. In addition, the bridge dynamic is expressed 
by using the second order differential equations in 
Eq. (9). In this study, the bridge dynamic is 
calculated by considering the first four-vibration 
mode and it is represented with four differential 
equations of second order. These four differential 
equations are transferred to eight equations by [16]. 
In order to solve this equation system, comprising 
of sixteen equations, a fourth order Runge-Kutta 
method has been used[16]. 

III. NUMERICAL ANALYSIS 

In this section, using four different VBA 
models (Model A, B, C, D) the effects of the 
absorbers on the vehicle and bridge dynamics have 
been studied extensively. System equations of 
motion Eq. (8a-d) for the VBA is solved by using 
Runge-Kutta algorithm of the fourth order, with a 
special m.file prepared MATLAB © environment. 
The integration time step size for numerical 
analysis during the study At = 0.01 s, the solution 
final time is taken to be t = L / v. In addition, the 
parameters used in this study for the systems VBA 
are presented in Table 1. 

When the vehicle moves over the bridge, 
it forces the bridge to vibrate. At certain vehicle 
speeds, the bridge enters into resonance and the 
amount of oscillation increases significantly. The 
velocity of the vehicle that causes resonance is 
called the critical velocity v ci and it is calculated 
using the circular natural frequencies of the beam. 
For a simply supported beam the natural 
frequencies of j th vibration mode is expressed 
as[17]: 

2 jYei 

w)= - —(rad Is) (11) 

mL 

When Eq. (11) and the excitation 
frequency of a moving vehicle on the bridge are 
rearranged, the following is obtained[16]: 

. .J/2 

w w vL&mo v 

a = —= -= —cA_± = —, (12) 

Wj 2 pfj j-p%EIa v cr 

Using Eq. (12), Fig. 2 shows the first three 
modes of vibration of the bridge beam. Frequency 
values of the first three modes of vibration of the 
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bridge beam are, respectively, /i=0.2108 Hz, 
/2=0.8432 Hz, and / 3 =1.8972 Hz. Using these 
values and Eq. (22), the critical speeds for vehicles 
moving on the bridge are calculated, respectively, 
as v cr i=42.15 m/s, v cr 2=168.62 m/s, and v cr3 =252.82 
m/s. 



Figure2.Mode shapes of the bridge. 


In order to analyze the effects of the 
absorbers on the dynamics of vehicles and bridges 
four different models of the VBA shown in Fig. 3. 
In the model A as shown in Fig. 3 a, an undamped 
mass spring system is suspended from the bridge at 
the mid-point. The main vibration displacement of 
the bridge system is affected mostly from the 
fundemantal mode of the beam, for such reason the 
absorber is placed at the mid point where the 
maximum deflection occurs at the midpoint when 
the beam vibrates with only the first mode. The 


first natural frequency of the bridge is// = 0.2108 
Hz. Therefore, the absorber is tuned to resonate at 
the first natural frequency 

/ = (2 p) 1 yjk j /m, = 0.2108 . Thus, it is expected 

that this absorber can reduce the vibration energy 
of the beam caused from the first vibration mode of 
the bridge. In the second model, model B, a 
viscous damper added to the spring mass system. 
The damping coefficient of the absorber was 
chosen as C2=10000 Nsm' 1 , not on purpose but to 
show the effect of any damping in the absorber. 
The critical damping value 

c cr2 = 2 m n w n = 2 s jk 2 m 2 = 26476.46 Nsm' 1 , and 
the damping ratio of the damper is 
z = c 2 / c cr2 = c 2 /(2m 2 w n )= c 2 f (2 ^k 2 m 2 ) = 0.377 
obtained using the parameters listed in Table 1. 
When this damper is added, the natural frequency 
of the absorber has reduced from 0.2108 Hz to 
0.195 Hz that is the damped natural frequency, 

where w d = w n -Jl- z 2 = 1.226 rad /s= 0.195 Hz. 
it is expected that adding damping may not be 
useful in order to reduce the vibration energy of the 
first mode. However, it should be known that in 
some case the damping could be considered in 
terms of design criteria. These issues are examined 
in detail in the following sections. 




Figure 3. Model of a vehicle-bridge-absorber interaction system (a) an undamped absorber at middle 
of the bridge; (b) A damped absorber at middle of the bridge; (c) three undamped absorbers at 0.25F, 0.5F, and 
0.75F of bridge length; (d) without absorber. 

The Models A and B can reduce oscillations originating from the first mode of the 
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bridge during vibration. Because these absorbers 
prepared at a middle point of maximum amplitude 
of the first mode of the bridge and their resonance 
frequency is set to the first oscillation frequency of 
the bridge. However, structures are not only 
affected from the first mode of vibration the 
contribution of the higher modes are also involved. 
Therefore, in addition to the first vibration mode 
for the second vibration mode of the bridge a 
second and a third absorber are placed at the 
locations 0.25 L and 0.75 L of the bridge; and they 
are tuned to the second natural frequency of the 
bridge beam. As seen from Figure 3 the peaks of 
the second mode shape of the beam are at these 
locations. In order to compare with the effects of 
the three previous models, the Model D has no 
absorber. The parameters of the bridge and the 
vehicle used in the study are presented in Table 1. 

Table l.Prpoperties of the vehicle and bridge. 


Bridge Vehicle parameters 


L( m) 

100 

m (kg) 

2172 

c 2 (Ns/m) 

10000 

E 

207 

* (kg) 

85439.6 

m/( Ns/m) 

10000 

(Gpa) 

0.174 

c (kg) 

2219.6 

m 2 ( Ns/m) 

10000 

/(m 4 ) 

20000 

kj (N/m 

280401. 

mj( Ns/m) 

10000 

E 

1750 

) 

59 

k 3 (N/m ) 

280401. 

(kg/m) 


&2 (N/m 

17525.0 


59 

Ceq 

(Ns/m) 


> 

8 




Based on the four different VBA models 
in Figure 3, Figure 4, displays the midpoint 
displacements Wb ( x=U2, t) of the bridge beam for 
a vehicle speed of v=90 km/h = 25 m/s that this 


x 10 3 



speed is chosen so that it is applicable in 
transportation. 

For Model D model that there is no 
absorber, the maximum midpoint displacement is 
20.34 mm when the vehicle is at about 0.75 L, 
while for Model A, the maximum displacement is 
observed as 15.8 mm, which means a decrease of 
22.32 percent when compared with no absorber 
case. After a damped absorber in Model B is used, 
the maximum midpoint displacement is 17.6 mm 
when the vehicle has reached to 81.2 percent of the 
bridge length. When a damping added to the 
absorber, the occurrence time of the maximum 
bridge midpoint displacement shifted to forward, 
because of the added phase difference from 
damping. Damping in the passive absorbers could 
not reduce the displacements much when compared 
to the undamped absorbers. For Model C, a special 
case when three absorbers are used, the maximum 
displacement of the midpoint of the bridge is 15 
mm and; that means the displacement has been 
reduced by 26.25%. All the results are presented in 
Table 2. 

Table 2.Comparation of the four different models according to 
bridge midpoint maximum displacement (m). 

Model Max. Mid-point Location of According to 

disp.(mm) the vehicle the Model D, 

(%) the reduction 


rate (%) 


A 

15.8 

73 

22.32 

B 

17.6 

81.2 

13.47 

C 

15 

87.5 

26.25 

D 

20.34 

75 

— 



Figure 4. The effect of the absorber upon bridge midpoint dynamic response for vehicle constant velocity v=25 
m/s and time step size At=0.01 s (a) Bridge midpoint displacement; (b) Absorber displacement. 


In traditional highway or railroad bridge 
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design, the Dynamic Amplification Factor (DAF) is 
considered. DAF means the ratio of the maximum 
dynamic displacement to the static displacement 
when the load is located at the midpoint of the 
bridge; and can be expressed as given below: 

,13) 

The DAF of any bridge system can be 
affected by many parameters such as vehicle 
velocity, mass of the vehicle, natural frequency of 
the bridge, and road roughness, etc. considering the 
Models A, B, C and D, the DAFs of the different 
velocities of the vhecile are presented in Figure 5. 
As shown in Figure 5, for the Model D, that is no 
absorber on the bridge the maximum value of the 
DAF is 1.65 for the vehicle speed v = 26 m/s, while 
for the Model B, the maximum value of the DAF is 
1.47 for the vehicle speed v = 18.5 m / s; and for 
the Model C the maximum value of the DAF is 
1.45 for the vehicle speed v = 13.5 m / s; and for 
the model A the DAF=1.37 at v=12.35 m/s. Up to 
a speed of 31 m/s, the Model A showed the best 
performance, but after this speed, the Model C was 
the best in reducing the vibrations. Considering 
this situation it can genarally be accepted that 
velocity of the vehicle determine the behaviour of 
the absorbers. This is because; the forcing 
frequency of the VBI system is determined by the 
velocity of the vehicle. 


IV. Conclusion 

In this study, in order to reduce the 
vibration of the vehicle that is caused by vehicle - 
bridge- interaction (VBI), the effects of various 
dynamic vibration absorbers on the bridge 
dynamics have been investigated. Due to the 
complexity of the loading in moving load problem, 
that the interaction is a non-linear coupled system 
of the two multiple degree of freedom subsystems. 
In the vibration of any beam system one can 
consider the general forced vibration response of 
the beam can be calculated by super positioning the 
mode shapes of consecutive modes at least 
including the first-four modes. However, in this 
assumption the effect of the first fundamental mode 
and the second mode is compromises the 80 
percent of the response even in very flexible 
systems. In such a case, the anti-node points of the 
two first nodes can be used as the location of the 
absorbers. The anti-node is the midpoint in the 
first mode, and they are at 0.25L and 0.75 L in the 
second mode shape. From the results of the 
analysis, the Model C in this study, represents this 
situation where three absorbers have been placed 
on the beam at the location of the anti-nodes of the 
first and second modes, is the best design for the 
reduction of the vibrations of the bridge. 
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